Landau-gauge condensates from the quark propagator on the lattice 
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We compute the dimension-2 condensate, {A^), and dimension-4 mixed condensate, (q^q), from 
the recent quenched lattice resuhs for the quark propagator in the Landau gauge. We fit the lattice 
data to the Operator Product Expansion in the "fiducial" region 1.2 GeV < Q < 3 GeV. Our 
result for the dynamical gluon mass at the scale of 10 GeV'^ is rriA = 600 — 650 MeV, in agreement 
with independent determinations. For the mixed Landau gauge condensate of dimension-4 we get 
as{qg^q) = (—0.11 ± 0.03) GeV^. This value is an order of magnitude larger than the (G^) gluon 
condensate. 
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I. INTRODUCTION 

The question of how constituent quarks arise dynami- 
cally has always been one of the most intriguing problems 
of QCD. The issue has prompted perturbative and non- 
perturbative approaches both in the continuum as well 
as on the lattice. Politzer H was the first one to com- 
pute the quark mass function using the Operator Prod- 
uct Expansion (OPE) in the high momentum regime in 
terms of the quark condensate (qq) in the Landau gauge. 
This calculation was corrected and extended to a gen- 
eral Lorentz gauge by Pascual and de Rafael The 
gauge-independent gluon condensate (G^) was included 
by Lavelle and Schaden where it was also forseen 
that a dimension-2 condensate, {A"^), should be present. 
Originally it was interpreted as a signature of sponta- 
neous gauge symmetry breaking. The dimension-4 mixed 
quark-gluon condensate, {q^q), was included in the anal- 
ysis of Lavelle and Oleszczuk ^] . 

More recently, Schwiiiger-Dyson approaches (for re- 
views see, e.g., Refs. |HIE0, 01 and references therein) 
have been intensely applied in an attempt to under- 
stand the nonperturbative physics in the infrared do- 
main. The phenomenological success of this approach has 
triggered a lot of activity on the lattice where the quark 
propagator has recently been computed after gauge fix- 
ing "Ciol m 113. However, the discussion of Ref. [H 
regarding the matching to the OPE is limited to the mass 
function. Remarkably, the dimension-2 {A-) condensate 
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is related to the quark wave function renormalization |3| . 

Early implications of a non- vanishing dimension-2 con- 
densate have been explored by Celenza and Shakin 
13j JJj|. More recently Chetyrkin, Narison, and Za- 
kharov [isl IT^ found that the inclusion of a tachyonic 
gluon mass parameter rriA ~ TOOMeV improves the phe- 
nomenology of the QCD sum rules in mesonic channels. 
For heavy quarks is proportional to the string tension 
of a short string, so it provides the short-range behavior 
of confining forces. Other phenomenological determina- 
tions of a non- vanishing gluon mass can be traced from 
the review (17i |. 

Although the dimension-2 {A"^) condensate naively 
breaks gauge invariance in the perturbative sense, a de- 
tailed analysis reveals that this is not so. As suggested in 
Refs. [iMll3il23i there exists a nonlocal gauge invariant 
condensate, 

which reduces to the (A^) condensate in the Landau 
gauge. Here g denotes the group element. A physical 
meaning has also been attached to this condensate by 
a perturbative gauge-covariant redefinition of the gluon 
field Further mounting evidence for the existence 

and physical relevance of the dimension-2 condensate 
in QCD has been also gathered from the lattice cal- 
culations [23 J analytic estimates jl^l, purely theoreti- 
cal considerations "M], and microscopic approaches [2^. 
Anomalous dimensions for the A^ condensate were cal- 
culated in Refs. [H HI [H IH H • 

The comparison of numerical lattice QCD calculations 
with analytic continuum approaches, such as the pertur- 
bation theory, the operator product expansion, or the 
Dyson-Schwinger approaches, requires a local gauge fix- 
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ing condition on the lattice. Thereafter it is possible to 
give a meaning of quark and gluon correlation functions. 
However, it is well known that there is no known local 
gauge fixing condition free of the Gribov copies (see, e.g., 
Ref. |32iJ and references therein). Therefore, one should 
keep in mind that when fixing the gauge there may still 
be differences in physical observables which become non- 
analytic functions of the coupling constant due to the in- 
fluence of the Gribov copies. If one restricts, however, to 
the fundamental modular region by a partial local gauge 
fixing, there may still be gauge-invariant operators under 
the residual subgroup and the BRST transformation [2^ . 
In the Landau gauge the only dimension-2 operator sat- 
isfying the above condition is precisely A^. 

In the present work we extract the dimension-2 {A^) 
condensate by comparing the lattice results for the 
quark propagator in the Landau gauge, presented in 
Ref. [13, m, to the OPE results of Refs. 0, ||. Our 
determination yields a novel estimate of the gluon mass, 
niA, as well as provides the first determination of the 
mixed dimension-4 condensate, (qg^q) (hereafter g de- 
notes the strong coupling constant). 
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II. LATTICE DATA FOR THE QUARK 
PROPAGATOR IN THE LANDAU GAUGE 



The inverse quark propagator can be parameterized as 
S~^{p) = ^A{p) — B{p), where A and B, dependent on 
the quark momentum, have the meaning of the vector 
and scalar quark self-energies. An equivalent parameter- 
ization is via the wave-function renormalization, Z, and 
the mass function, M, defined as 



S{p) 



Zip) 



i) - M{p) ' 



Z{p) - A-\p), M{p) = B{p)IA{p). 



(2) 
(3) 



The quark propagator was calculated in Landau gauge 
using the "Asqtad" improved staggered action. The 
gauge ensemble is made of 100 quenched, 16'^ x 32 lattices 
with a nominal lattice spacing a — 0.124 fm, set from the 
static quark potential. This data was first published in 
Ref. [13. 

The results for M and Z as functions of the Euclidean 
momentum Q are shown in Fig. ^ at various values of 
the current quark mass m. The data for M asymptote 
at large Q to the value of m, indicated by the horizontal 
lines. We note that the data at highest values of Q are 
not perfect, with some visible wiggles and a tendency of 
falling off at the end, which may be attributed to the 
finite-size effect. Yet, up to Q ~ 3 GeV the tails in 
M and Z look very reasonably, reaching plateaus before 
"hitting the wall" . 

As a matter of fact, the tail in M in the "fiducial" 
region of 1.9 GeV < Q < 2.9 GeV was used successfully 



FIG. 1: The quark mass function M (top) and the wave- 
function renormalization Z (bottom), plotted as functions of 
the Euclidean momentum Q. The data comes from quenched 
lattice calculations in the Landau gauge of Ref. i. Various 
sets of points correspond to the current quark masses m = 
29, 42, 54, 80, 105, 150, 225, and 295 MeV, indicated by 
horizontal lines in the top panel. In both panels the highest 
sets of points correspond to highest values of m. 



in Ref. 12] to verify the expression 



M(Q) 
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where oJm = 12/(33 — 2A^/) with Nf = flavors, {qq)fj. 
and m(fi) are the quark condensate and the current quark 
mass at the scale fi, respectively, and Aqcd — 691 MeV 
in the MOM scheme. This shows that the data is accu- 
rate enough to be verified against the perturbative QCD 
predictions. 

The data for Z{Q) from Ref. are shown in the bot- 
tom panel of Fig. ^ A very weak dependence on m has 
been noted, except perhaps at low Q. Asymptotically, 
Z{Q) — > 1, as requested by the canonical normalization 
of the quark fields. At lower values of Q the departure 
of Z from unity is sizeable, with a long-range tail clearly 
visible. 



III. MATCHING OPE TO LATTICE DATA 

In our further analysis we will work with the function 
A{Q). The data for the vector quark self-energy A{Q) 
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FIG. 2: The fit to the tail of the vector quark self-energy, 
A{Q). The solid line corresponds to the asymptotic formula 
Q with the optimum parameters (|SJ, while the data (in- 
cluding all values of the current quark masses m) are taken 
from the quenched lattice calculation in the Landau gauge of 
Ref. 0. The asymptotic curve is drawn in the fiducial region 
of 1.2 GeV < Q < 3 < GeV. 
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FIG. 3: The contours corresponding to the fit of the data 
of Ref. i for the vector self-energy A{Q) in the Landau gauge 
to formula l^l. The dot indicates the optimum values of 
Eq. ©. The curves correspond to 68% and 95% confidence 
levels. Note a large correlation between C2 and a. 



may be parameterized at sufficiently large values of Q as 



A{Q) 



1 



C2 



C4_ 
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In the fitting procedure we must decide on the matching 
region in Q. Certainly, this choice will affect the results, 
yielding a systematic error. The values of Q cannot be 
too large due to finite-size effects, nor too small, where 
the expansion ^ is no longer accurate. We perform the 
fit in the range 1.2 GeV < Q < 3 GeV, which yields 
the optimum values 

C2 = (0.37 ±0.04) GeV^ C4 



-0.25 ±0.06) GeV^ (6) 



The errors have been calculated by jackknife. The value 
of x^/DOF is 0.51, but one can see from the plot of Fig. 13 
that there is a sizeable correlation between C2 and C4. 

We have also performed the fit with three terms, ap- 
pending formula with the term cg/Q^. The result is 



C2 = 0.39 GeV^ a = -0.37 GeV*, and cq = 0.15 GeV^ 
with errors that overlap with the central values for C2 and 
C4 of Eq. However, due to large correlations between 
C4 and cg, no reliable information may be extracted from 
this three-parameter fit. More accurate data and a larger 
range of momenta will allow for a better determination 
of the l/Q^-expansion. 

Next, we will compare the obtained values of Eq. © 
to theoretical predictions and extract estimates for the 
Landau-gauge condensates. At D = A the vector self- 
energy read out from the propagator of Ref. 0, Q is 



A{Q) 



1 



(7) 
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where fj, denotes the renormalization scale, 
to Eq. ® we find for three colors 

asin^){A^)^ = (0.36 ±0.04) GeV^, 



Comparing 



(8) 



or 



and 



g'^(A^) = (2.1 ±0.1 GeV)^ 



(9) 
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(-0.11 ±0.03) GeV". 



Since (S^G^) 



(10) 

0.009 GeV" m, the contribu- 
tion of the (G^) condensate to Eq. pO|) is negligi- 
ble compared to the mixed-condensate term. Thus 
a,(At2)(gg4g)^ = (-0.11 ±0.03) GeV^ 

The errors quoted throughout the paper are statisti- 
cal. In addition, there are certain systematic errors orig- 
inating from the choice of the fitted function A{Q) of 
Eq. and from the choice of the fiducial region in Q. 
Quantities quoted in physical units are also subject to 
the uncertainty in scale that goes with quenched lattice 
simulations. 



IV. 



COMPARISON OF (A^) TO OTHER 
DETERMINATIONS 



The Landau-gauge condensates considered in this pa- 
per are not renorm-invariant quantities, thus their values 
evolve perturbatively with the scale. The QCD evolution 
for (^2) has been worked out in Ref. [H ll IH H |3ll , 
with the leading-order result 



as{^i'){A')^^a,i^i'y-'-^^f'\ 



(11) 



where 7^2 =35/4 and /3o = 11 correspond to evolution 
with no flavor. We use agip"^) = 47r/(9 log[^^/A^]), with 
A = 226 MeV for the evolution at the leading order. The 
exponent in Eq. Hll|l is equal to 9/44, hence the evolution 
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is very slow. For instance, the change of //^ from 1 GeV^ 
up to 10 GeV^ results in a reduction of as(A^) by 10% 
only. 

Most estimates in the literature refer to the gluon mass, 
related to the (A^) by the formula m\ = ^g^(A^). 
Our estimate ((HJ, when evolved with Eq. from 2 
to 10 GeV^, yields 

ruA = (625 ± 33) MeV. (12) 

Evolution from 1 to 10 GeV^ gives tua = (611±32) MeV, 
while evolution from 4 to 10 GeV^ produces ruA = (635 ± 
34) MeV. These values are close to many estimates made 
in other approaches. In particular, most of the numbers 
listed in Table 15 of the review of Ref. ^3 obtained 
by very different techniques are in the range 0.5-1.5 GeV. 

V. CONCLUSIONS 

We have attempted to match the OPE to the quenched 
lattice data for the vector quark energy in the Landau 



gauge. The obtained value of the dimension-2 Landau- 
gauge gluon condensate, {A"^), of Eq. |(SJ) and the cor- 
responding estimate for the gluon mass of Eq. (|12|l are 
consistent with other estimates in the literature. Thus 
the lattice provides an independent way of determining 
this condensate. The estimate for the dimension-4 mixed 
quark-gluon condensate of Eq. (|10(l , made to our knowl- 
edge for the first time, is an order of magnitude larger 
compared to the (G^) condensate. 
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